Two-body transients in coupled atomic-molecular BECs 
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We discuss the dynamics of an atomic Bose-Einstein condensate when pairs of atoms are converted into 
molecules by single-color photoassociation. Three main regimes are found and it is shown that they can be 
understood on the basis of time-dependent two-body theory. In particular, the so-called rogue dissociation 
regime [Phys. Rev. Lett., 88, 090403 (2002)], which has a density-dependent limit on the photoassociation 
rate, is identified with a transient regime of the two-atom dynamics exhibiting universal properties. Finally, we 
illustrate how these regimes could be explored by photoassociating condensates of alkaline-earth atoms. 



The conversion of atom pairs into molecules, using ei- 
ther Feshbach resonances jl| or photoassociation fl, serves 
as a tool to probe the many-body properties of ultracold 
gases l^. In particular, photoassociation, the process of as- 
sociating atoms with a resonant laser light, was recently used 
to observe pair correlation in a one-dimensional Bose gas J3l 
and the crossover between Bose-Einstein condensate (BEC) 
and Bardeen-Cooper-Schrieffer (BCS) superfluid [5]. Con- 
versely, it can be used to reach new regimes. Many-body 
theories have suggested the coherent conversion of an atomic 
BEC into one of molecules macroscopic superposi- 

tion (Hi, and production of correlated atom pairs at high laser 
intensity J3, 3 HI- Several experiments have made the first 
steps in these directions il2lll3lll4ll5ll , but have been lim- 
ited by inherent losses or insufficient laser power. Refer- 
ence 111 111 identified three regimes of photoassociation as a 
function of loss and laser intensity. The intriguing density 
dependence of the regime boundaries suggested that they are 
associated to many-body effects. 

In this letter, we first apply time-dependent two-body the- 
ory to photoassociation with a single continuous laser and 
distinguish three transient regimes. We then show how these 
transients explain the previously identified regimes in the 
many-body theories. 

For two atoms of mass M interacting with a resonant laser, 
photoassociation is described by two equations coupling a 
scattering and a molecular channel {lfjft . 



ihcf>(r)= - — V 2 + U(r)j <j>{¥) + W(r)cb m (r) 

iH m {r)=( K - ^V 2 + U m {r) - i^)<p m (f) + W(r)<j>(r), 

where r is the relative separation of the two atoms, 4>(r) and 
4> m (r) are the components of the relative motion wave func- 
tion for the scattering and molecular channels, U and U m are 
the interaction potentials in each channel, 7 is the sponta- 
neous emission rate from the molecular channel (we assume 
that decayed molecules are lost from the system), and W 
couples the two channels. W is proportional to the square 
root of the laser intensity. We expand \<p) and \(j> m ) in the 
bases of eigenstates of — fjV 2 + U and -fjV 2 + U m , re- 



spectively, and assume that only the scattering eigenstates 
|<p^) (indexed by wave vector k) are relevant in the scatter- 
ing channel, and that a single bound eigenstate | </?,„) is reso- 
nant in the molecular channel. Choosing (<f m Wm) = 1 and 

(^fcl^p) = (2n) 3 5 3 (k — p), one obtains 



ihCAt)=E k CM) +wtC m {t) 



(1) 



ihC m {t) = (A - n/2)C m (t) + iTT^w^t), (2) 
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where and C rn are the amplitudes in states \(p?) and \tp m ), 
Ek = h 2 k 2 /M, A is the resonant bound state energy with 
respect to the scattering threshold (which can be adjusted by 
tuning the laser frequency), and wj: = (ip m \W\(pj:} are the 
coupling matrix elements. According to Wigner's threshold 
laws, wj: goes to a constant w for low k <C 1 / 1 r c | , where r c 
is the largest of the extent of the molecule, the van der Waals 
length [2], or the scattering length a associated with U. 

In ultracold gases, atoms collide at nearly zero energy. For 
the stationary solution at zero energy, Cg goes to — AnA/k 2 
for low fc, where A is the optically-induced complex scatter- 
ing length tHinl given by 4-Kh 2 A = -M M 2 /(A - A' - 

The imaginary 
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17/2) and the light shift A' = / ^ 
part of A is related to the loss rate coefficient 



Si:h J . 2 T 
K = — —ImA = -Im— 
M h A 



i 7 /2' 



(3) 



which corresponds to the number of atoms lost per unit of 
time and volume due to photoassociation to the excited state 
and subsequent decay by spontaneous emission. From ki- 
netic theory the density p of remaining atoms in a thermal 
gas is expected to follow the rate equation 



-Kp 2 . 



(4) 



We now take into account the fact that the laser is turned 
on at t = 0, creating a strong perturbation to the two-body 
system. As a result, transient regimes appear before the sta- 
tionary solution is reached. Initially, the two atoms are in 
the scattering channel with nearly zero collision energy, i.e. 
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Cj: = (27r) 3 <5 3 (fc) and C m = at t = 0. For i > Owe 
choose to decompose CAt) as follows 

Cg(t) = (2^) 3 <5 3 (fc) + Cf (t) + Cf "(t), (5) 

where C^ d (t) = —-^■C m {t) is the adiabatic response to the 
turn-on of the laser (obtained by setting G? = in Eq. (0), 
and (7^ yn (t) is the dynamic response. Solving for ci vn from 
Eqs. (0 and © and inserting it into Eq. 0, we obtain 



zftC m (i) = (A- A'-il)c m (t) 



1 - i 
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where we have set wj: = w, which is valid for t t c , t c = 
Mr\ /h being in most cases a very short time scale on the or- 
der of 10 ns. We then assume that C m (t) is localized at short 
times, which leads to the Ansatz f Q d,TC m (r) / y/t — r = 
oiC m (t) I y/t, where a is to be determined, and we approxi- 
mate ihC m (t) by the short-time expression ihC m (t)/t. This 
leads to 



C m {t) 



A - A' - 



(6) 



where S(t) = 
a* (l 



and f 



S{t) ~ f 

are time-dependent 
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shifts and broadenings of the molecular level, and t w = 
i^ 3 (^) 4 . In the limit of large S(t), one has C m (t) = 
~ V^wt which justifies that C m (t) oc t~ x / 2 is local- 
ized at short times and yields a = ir/2. For small S(i), our 
Ansatz may not hold but this has little consequence precisely 
because S(t) is small. 

We can now calculate an instantaneous rate coeffi- 
cient, based on the time variation of the population 

Ik~o (t J)* l^fcl 2 ^ n ^ e initial state. This population is always 
inifinite, because we started from a state which is not normal- 
izable over an infinite volume. However, from Eq. (Q3, its 
time derivative has a finite value, which we identify as minus 
the instantaneous rate coefficient, —K{t). Using Eq. dD, we 
can simplify it to 



K{t) 



-lm(w*C m (t)). 



(7) 



Depending on the relative strength of the terms in the de- 
nominator of Eq. ©, this coefficient goes through three sub- 
sequent regimes illustrated in Fig. [T] linear (a), square root 
(b), and constant with time (c). Namely, 



for t c < t < t. 



(8) 



(b) K(t) = 4/tt (h/M) 3/2 Vt, foit w < t < t A (9) 



(c) K(t) = -Im 



h A-A'-*7/2' 



fort A <t. (10) 
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Figure 1: Schematic evolution of the rate coefficient (thick black 
curve) showing the three regimes (a) (b), and (c) at sufficiently 
high laser intensity (t w i>0- The grey curves correspond to the 
limiting expressions J81 1 Ob . In particular, the parabolic curve asso- 
ciated with Eq. ((9|( is a universal upper limit to the rate coefficient. 



where t A = a 2 ^-|^4| 2 . Note that t w < t A only for 
high laser intensity. If t w > t A , regime (a) occurs for 
J < f 7 and regime (c) occurs for t ^> t 7 , where t 7 = 

^/v/(A-A') 2 + ( 7 /2) 2 . 

Photoassociation converts the initial state into the molec- 
ular state and the loss rate coefficient grows linearly in 
regime (a). In regime (b) the laser also drives molecules back 
to the atom-pair continuum. This broadens the resonance 
and the loss rate coefficient still increases but more slowly. 
Finally, in regime (c) the molecules start decaying sponta- 
neously, and the rate coefficient reaches its steady-state max- 
imal value, which is the rate coefficient Eq. (0 obtained from 
the stationary solution. 

The atomic density is expected to follow Eq. (0]i with the 
time-dependent rate coefficient (0. Whether the transient 
regimes (a) and (b) are observable depends on whether an 
appreciable fraction of the density is photoassociated over 
these time scales. Let us define the depletion time t p for 
which an appreciable fraction is depleted, i.e. K(t p )pot p m 
1, where po is the initial density. Three cases are possible. If 
t A , t-y t p , then only the final constant rate coefficient is 
relevant. If t w <C t p -C t A , then the first regime to lead to 
observable losses is regime (b). If t„ < t w , i 7 , then the only 
relevant regime is (a). 

Interestingly, when regime (b) dominates, the molecular 
amplitude vanishes and the system has a universal behavior. 
The loss rate (0 and the atom pair distribution 



CAt) = (2^) 3 <5 3 (fc)+ 
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do not depend on the microscopic details of the transition, 
but just on the mass of the species. Here, Erf denotes the 
error function. The condition t w <i p < t A , needed for the 
observation of this regime, is equivalent to 

pl'\A\ > (2tt)- 2 / 3 and p\j Z B « (2tt)- 2 / 3 (11) 

where the length B is 8 (^jj^j ■ Finite temperature adds 

the condition that K{t p ) be smaller than the unitarity limit 
jjA, where A is the de Broglie wavelength. It fortuitously 
coincides with the condensation condition for bosons pA 3 > 

We now turn to a many-body description. Photoassocia- 
tion in a uniform BEC can be described (up to first order in 



3 



■o 

c 

cn 
c 

'c 
Q 



(c) 
adiabatic 
regime 



/ (a) 
coherent 
regime 



(b) 
rogue 
dissociation 
regime 



Laser intensity 7 

Figure 2: Regimes of photoassociation in a BEC as a function of the 
"detuning and width" ^/(A — A') 2 + (7/2) 2 and the laser inten- 
sity I, for a fixed density. The boundaries are indicated by dashed 
lines and equalities of time scales defined in the text. 



a cumulant expansion 118]) by the three equations iflolfTlll 
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» = r J \ U(r)$(r) + W(r,t)$ m (r 
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U(r))$(r) + W(r)$ m (r) + 2im^ 
U m {r)- i 7 /2W„(r) + W(rMr) 



where ^ is the condensate wavefunction (here a complex 
number), and $(f) and $ m (f) are the pair wavefunctions 
in the scattering and molecular channels. Higher-order cu- 
mulants (such as the normal density of noncondensate atoms 
lisin contribute significantly only for t > t p , and do not 
change the dynamics up to 10 /is in the examples to follow. 
Inelastic collisions B2011 . not included here, do not affect the 
atomic BEC during this time frame for typical rate coeffi- 
cients ~ 10~ 10 cm 3 s _1 . 

As in the two-body case, we can write 



$(f) = 



■(CI 



ad 
h 



s-idyn 



where the adiabatic part C2 



ad 



i-K-* m +5fef 2 ) and 



g? = (0\U\(p^), |0) being the zero-momentum plane wave. 
Elimination of the adiabatic part is now crucial, because it 
introduces the light shift A' but also the coupling constants 
u>£ and g? without resorting to contact interactions 111 ill . We 
find 



(2tt) 3 



(12) 



d 3 k 



dyn 



13) 



ihct yn = E k ci yn - ihdf. 

k k k 



(14) 



These equations are similar to those of Ref. [9], but contain 
no ultraviolet divergence. Thus, we can thus safely set = 
w and gj: = g = 4irh 2 a/M without any renormalization. 
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Figure 3: (color online) On-resonance photoassociation (A — A' = 
0) of a condensate of sodium (a), ytterbium (b) and strontium (c) 
for typical transitions. In each row, the right panel shows the pop- 
ulation evolutions, according to Eqs. d 1 21 14b . based on the pa- 
rameters indicated by the black dot in the left panel, which is a 
regime diagram similar to Fig. [2] Solid curve: atomic conden- 
sate population; dot-dashed curve: molecular population; dotted 
curve: correlated pair population. The short-dashed curve shows 
the atomic condensate population following from Eq. (O with 
the rate J 1 6b - For sodium, we used the conditions of Ref. fl2h 
(|ImA|/J = 2.95 fm/(WcirT 2 ), j/H = 18 MHz). Typ- 
ical intercombination transition parameters were used for ytter- 
bium |H QlmAUI = 2.12 nm/(Wcm~ 2 ) , j/h = 364 kHz) 
and strontium lHj (\lmA\/I = 2.12 nm/(Wcm" 2 ), j/H = 15 
kHz). Thus, we have {t w , t y , t A ) = (0.17, 0.017, 0.0036) /is for 
sodium, (0.044, 0.88, 35) ps for ytterbium, and (116, 21, 7.8) ps 
for strontium. For all cases, the initial density is po = 6 • 



10 14 cm^, andt D 



5 ps. 



When the dynamical part C k yn is negligible we obtain the 
familiar set of coupled Gross-Pi taevskii equations introduced 
in Ref. 112111 . If we assume that = ^/po and ^ m = ini- 
tially, these equations admit two limiting regimes ifTH |2lll . 
In the adiabatic regime, w 2 p < (A - A') 2 + (7/2) 2 , 
the molecular wave function can be adiabatically eliminated, 
and the condensate density p = |^| 2 then follows Eq. (O 
with the rate coefficient ([3]l predicted by the stationary two- 
body theory. On the other hand, in the coherent regime, 
w 2 pa > (A - A') 2 + (7/2) 2 , vf/ and * m exhibit coherent 
Rabi oscillations at a frequency £1 = w^fpo/h. 



The case where the dynamical part C- 



dyn 



cannot be ne- 



glected corresponds to the "rogue dissociation limit" of 
Ref. [9]. Ref. Il 111 showed that this occurs when 

m, . 1 /TTX 1 / 3 1/3 1 /7T\ l/ 3 

py\ A \~>2\2) and Po L> 2\2) (15) 
for the adiabatic and coherent regimes, respectively, where 
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was identified with a many-body length. Fig- 



ure|2]shows the three regimes: adiabatic, coherent and rogue 
dissociation. Most intriguing has been the dependence of the 
regime boundaries on the density. In particular, increasing 
the density makes it easier to reach rogue dissocation from 
the adiabatic regime, but more difficult from the coherent 
regime, according to Eqs. < fT3T >. 

In light of the previous two-body analysis, we can now 
interpret the adiabatic, coherent and rogue regimes in terms 
of the two-body regimes (a), (b) and (c). Using the same 
approximations as in the two-body theory, we can reduce the 
many-body equations (1121114b to a rate equation with a time- 
dependent coefficient 



K(t) 



2 Im M 2 +g-S(*) 



S(t) 



(16) 



We emphasize that this rate coefficient does not depend 
on the density and originates essentially from the time- 
dependent two-body coefficient (0. It follows that the three 
regimes (a), (b) and (c) also apply for ( [TBI ). In fact, the con- 
ditions for the observation of regime (a), (b) or (c) during 
the depletion time t p , expressed in terms of t w , t~ ( , and tA 
are equivalent to the boundary conditions for respectively the 
coherent, rogue dissociation, and adiabatic regime, as shown 
in Fig. |2 Thus, the density dependence of these boundaries 
originates simply from the usual rate equation of kinetic the- 
ory. Note in particular that the rogue regime boundaries ( fT5l ) 
are equivalent to the conditions ( fTTT i for the observation of 
the universal regime (b). This shows that the length L has in 
fact no special significance. 

The right panels of Fig. [3]compare the atomic condensate 
evolution from Eqs. ( 1121114b with that following from the rate 
coefficient Eq. ( TTol ) in the three regimes. It shows that the 
short-time evolution is always consistent with two-body dy- 
namics and kinetic theory, justifying the approximations we 
have made. Only in the coherent regime, for t > t p , is the 
condensate nature of the gas revealed due to collective effects 
which cannot be described by a rate equation. 

For alkali-metal atoms experimentally studied in 112, 151. 
the molecular state has a short lifetime on the order of 10 ns. 
As a result, these experiments have been confined to the adi- 
abatic regime. Figure^ shows a typical case for sodium. To 
reach other regimes, very large intensities are needed. On the 
other hand, photoassociation near narrow intercombination 
lines leads to much longer-lived molecules. Figures^ and 
[3J; show on-resonance photoassociation in condensates of yt- 
terbium and strontium, for typical states below the intercom- 



bination line. For moderate intensities, it appears possible to 
reach the universal regime (b) of pair dissociation. This cre- 
ates correlated pairs of atoms from a condensate, analagous 
to correlated photons ll24ll . For strontium, it should be pos- 
sible to reach the coherent regime, at least partially. Note, 
however, that inelastic collisions, while not significantly af- 
fecting the atomic condensate, reduce the molecular and non- 
condensate populations by about a half in Fig. [3};. 

In summary, we showed how different regimes of photoas- 
sociation in a BEC originate from transient shifts and broad- 
enings in the two-atom dynamics. These have simple ana- 
lytical expressions, and lead to a universal behavior in the 
rogue/pair dissociation regime. 
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